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ABSTRACT 



Context. Based on the analysis of a small sample of BATSE and Konus gamma-ray bursts (GRBs) with know redshift it has been reported that 
the width of the autocorrelation function (ACF) shows a remarkable bimodal distribution in the rest-frame of the source. However, the origin 
, of these two well-separated ACF classes remains unexplained. 
t-H ' Aims. We study the properties of the bursts belonging to each ACF class and look for significant differences between them. 
CD Methods. We complement previous ACF analysis studying the corresponding power density spectra (PDS). With the addition of Beppo-SAX 
t data and taken advantage of its broad-band capability, we not only increase the burst sample but we extend the analysis to X-ray energies. 

' Results. The rest-frame PDS analysis at y-ray energies shows that the two ACF classes are not simply characterised by a different low 
(-H , frequency cut-off, but they have a distinct variability as a whole in the studied frequency range. Both classes exhibit average PDS with 
O |' power-law behaviour at high frequencies (/' > 0. 1 Hz) but significantly different slopes, with index values close to those of Brownian (-2) 
and Kolmogorov (-5/3) spectra for the narrow and broad classes respectively. The latter spectrum presents an additional PDS component, 
a low-frequency noise excess with a sharp cut-off. At X-ray energies we find the power-law index unchanged for the broad class, but a 
significantly steeper slope in the narrow case (~ -3). We interpret this as an indication that the broad class bursts have weaker spectral 
^ evolution than the narrow ones, as suggested also by our analysis of the ACF energy dependence. The low and high frequency PDS components 
• • ■ may then arise from two radiating regions involving different emission mechanisms. We compare our GRB sample conditioned by afterglow 
^ ' detections with a complete, flux limited BATSE sample, finding a significant bias against narrow ACF bursts. 

Key words, gamma rays: bursts - gamma rays: observations - methods: data analysis - distance scale 
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1. Introduction 

Determining the relevant timescales for any astronomical phe- 
nomenon is essential to understand its underlying physical pro- 
cesses. However, in spite of extensive research, temporal stud- 
ies on the prompt emission phase of long gamma-ray bursts 
(GRBs) are not yet able to describe and explain their basic tem- 
poral properties. The main challenges encountered in the tem- 
poral analysis of GRBs are related to intrinsic characteristics 
of the emitted signal. Firstly, bursts are non-repetitive short- 
term events. Consequently, the total duration of the emission 
in a given observational energy window is the first timescale 
used to characterised them. Through out this paper we will 
only consider the class of long bursts, i.e., those with dura- 
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tion time Tgo > 2 s. Second, burst light curves (LCs) show a 
remarkable morphological diversity and they appear to have a 
composite structure. While a significant fraction of bright long 
bursts (~ 15%) exhibits a single smooth pulse structure, in most 
cases they appear to be the result of a complex, seemingly 
random distribution of several pulses. Burst pulses are com- 
monly described as having fast-rise exponential-decay (FRED) 
shape, although the decay is not strictly exponential. Therefore, 
the second timescale that seems relevant for the description 
of a burst is a "typical" pulse duration. However, analysis of 
the pulse parameters has shown broad log-normal distributions 
not only among different bu rsts, but also within a single burst 
(see, e.g jNorris et al.ll996h . Consideration of oth er timescales 
migh t be relevant, e.g., for GRBs with precursors (Koshut et al. 
1995) or when long quiescent periods occur (Nakar & Piran 



2 



L. Borgonovo et al.: Variability classes in GRBs 



20021) . although these temporal features appear only in a small 
fraction of bursts. 

Due to these characteristics, much of the GRB temporal 
analysis has been done directly on the LCs, i.e., modelling the 
pulses and studying their shape and distribution. Standard lin- 
ear analysis tools must be used with some caution, since most 
of the inferences based on them would in principle require 
"long" stationary signals, i.e., the most suitable bursts are those 
where the duration is much longer than the typical pulse width. 
This is the case of the temporal analysis based on power den- 
sity spectra (PDS). Individual PDS of GRB have very diverse 
shapes, and they do not seem to have common features, al- 
though the longest bursts show spectra with a consistent power- 
law behaviour. One way to overcome these limitations is to es- 
timate an average PDS from a sample GRBs. This approach 
will only produce physically meaningful results if each burst 
can be considered a realisation of the same stochastic process, 
i.e., there are no subclasses in the sample. 

Under th is assumption Beloborodov. Stern. & Svenssonl 
dl998l l2000h calculated with a large sample of bright GRBs, 
an average PDS showing a clear power-law behaviour extended 
over two frequency decades (approximately within the 0.01- 
1 Hz frequency range), and more remarkably with an exponent 
value approximately equ al to that of the Kol mogorov spectrum 
found in fluid turbulence dKolmo goro vl 1 94 1 1) . The significance 
of a thus obtained average PDS depends on two important ad- 
ditional factors. First, light curves have to be normalised to bal- 
ance the weight between bursts of different brightness. It is not 
clear at this point which norm should be used to pr oduce the 
most meaningful average. iBeloborodov et al. I d 19981) favoured 
the use of the peak flux normalisation, however they tested 
several other norms with qualitatively similar results. Thus, 
a good normalisation should give better convergence but the 
norm should not affect the final result for a sufficiently large 
sample. The second problem is the shift in frequencies pro- 
duced by cosmic dilation eff ects. Having no redshifts z deter- 
mined for their GRB sample, Bel oborodov et alJ (2000) did not 
correct for these effects. However, they argued that if the under- 
lying PDS shape for every GRB is a featureless single power- 
law with a constant exponent, the frequency shifts will not af- 
fect the obtained average PDS. If this is the case, considering 
that the two decades range of the power-law is much larger than 
the standard deviation of the redshift distribution (o% ~ 2 based 
on the few known redshifts), indeed the shifts should just smear 
the cut-off frequencies. 



The same statistical approach was used by Fenimore et al 



( 1995, hereafter F95) in their study of the average autocorre- 
lation function (ACF) of a sample of bright GRBs. The ACF 
gives a measure of the correlation between different points in 
the light curve that are separated by a given time lag. Since 
it is the Fourier transform of the PDS, it contains in principle 
the same information that can be visualised in a different way. 
Therefore, the same caveats regarding the average PDS apply 
to the average AC F. It was only aft er the discovery of the af- 
terglow emission dCosta et alJI 1997D and the determination of 
their redshift for a significant number of bursts t hat w e were 
able to address some of those issues. iBorgonovol ( 2004 , here- 



after B04) showed for a sample of 16 bright GRB with known z 



that when corrected for cosmic dilation effects the ACFs exhibit 
a clear bimodal distribution. Using as a measure the half-width 
at half-maximum, there is a highly significant gap between a 
narrow and a broad width class, the separation in standard devi- 
ations being >7<x. The estimated local or intrinsic values (i.e., 
those calculated at the rest-frame of the source) for the average 
widths were 1.6 s and 7.5 s, and the relative dispersions were 
32% and 4% for the narrow and broad classes, respectively. It 
is remarkable the low dispersion found in the last subset, which 
comprised ~ 1/3 of the total sample. 

This article builds on the ACF analysis done in IboI . In 
Sect. [2] we present our data samples and in Sect. [3] we briefly 
describe the methods used in the subsequent temporal analysis. 
In Sect. H] we strengthen previous findings based on the ACF, 
expanding the previous sample of GRB with known z by the in- 
clusion of proprietary data from the BeppoSAX mission, and we 
complement the temporal analysis estimating the intrinsic PDS 
for each of the subsets identified using the ACF. Furthermore, 
using the broadband capability of BeppoSAX combined instru- 
ments, in Sect. [5] we are able to extend the study to the X-ray 
energies. In Sect. |6l we investigate whether the typical values 
of several physical parameters commonly used to characterised 
GRBs differ significantly between the two temporal classes. We 
look into the problem of the energy dependence of the ACF and 
discuss possible biases in our sample of bursts with known z in 
Sect. |7] In Sect. [8] we discuss our main results. 

2. Data 

This work is mainly based on the analysis of light curves 
from GRBs with known redshift. Given the scarce number of 
cases available for study, we combined data (in the gamma en- 
ergy band) from three instruments to improve our statistics. 
Increasing a sample in this way presents an obvious trade-off, 
since we use count time series and the difference between in- 
strument responses introduce an additional dispersion that may 
counteract the benefits. For this reason, we initially analysed 
the data of the bursts that were observed by more than one in- 
strument, evaluating whether the differences were acceptable 
for our purposes. 

The comparisons were made taking BATSE as the refer- 
ence instrument. Its data comprise half of our GRB sample 
(in the gamma energy band), showing the best signal-to-noise 
ratio (S/N) with relatively low directional dependence thanks 
to the large collecting area of its eight Large Area Detectors 
(LADs) placed on eac h corner of the Com pton Gamma-Ray 
Observatory (CGRO; Fishman et al. 19891) . giving full sky 
coverage. It flew during the period 1991-2000 collecting the 
largest GRB catalog up to date. The CGRO Science Support 
Center (GROSSC) provides the so-called concatenated 64 ms 
burst data, which is a concatenation of the three standard 
BATSE data types DISCLA, PREB, and DISCSC. All three 
data types have four energy channels (approximately 25-55, 
55-110, 110-320, and > 320 keV). The DISCLA data is a 
continuous stream of 1.024 s and the PREB data covers the 
2.048 s prior to the trigger time at 64 ms resolution, both types 
obtained from the 8 LADs. They have been scaled to overlap 
the DISCSC 64 ms burst data, that was gathered by the trig- 



L. Borgonovo et al.: Variability classes in GRBs 



3 



Table 1. Sample of 22 GRBs with known redshift. The columns give the name of the GRB, the source instrument for the y-ray 
data, the X-ray instrument when available, the measured redshift z, the corresponding reference, the ACF half-width at half- 
maximum in the y-ray energy band w y , the half-width corrected for time dilation w' y , the corresponding half-widths wx and w' x 
from the X-ray light curves when available, the width ratio between the two energy bands wx/w y , the index £ assuming a width 
energy dependence W y (E) cc E~^, and the ACF width class. 



GRB Instrument (y) 


Inst. (X) z 


Ref* 


Wy 


K 




w 'x 




$2ch 


Class rf 


(55-320 keV) a 


(2-28 keV) 




(*) 


(s) 


0) 


0) 









970228 


GRBM 


WFC 


0.695 


(1) 


1.3 ± 


0.1 


0.77 ± 0.06 


3.5 : 


t 0.2 


2.1 : 


±0.1 


2.7 






n 


970508 


BATSE/GRBM 


WFC 


0.835 


(2) 


2.7 ± 


0.1 


1.47 ± 


0.05 


10.3 


± 0.4 


5.6: 


±0.2 


3.8 


0.06 


±0.05 


n 


970828 


BATSE 




0.9578 


(3) 


15.33 i 


:0.06 


7.83 ± 


0.03 












0.09 


±0.03 


b 


971214 


BATSE/GRBM/Konus 


WFC 


3.418 


(4) 


8.02 ± 


0.08 


1.81 ± 


0.02 


11.8 


±0.8 


2.7 : 


±0.2 


1.5 


0.10 


±0.04 


n 


980326 


GRBM 


WFC 


1.2 


(5) 


1.34 i 


:0.1 


0.61 ± 


0.04 


2.35 


± 0.4 


1.05 


± 0.2 


1.7 






n 


980329 


BATSE/GRBM/Konus 


WFC 


3 ± 1 


(6) 


5.96 ± 


0.02 


1.5 ± 


0.5 


8.8 : 


±0.4 


2.2: 


±0.7 


1.5 


0.05 


±0.01 


n 


980425 


BATSE/GRBM 


WFC 


0.0085 


(7) 


7.62 ± 


0.08 


7.56 ± 


0.08 


14.8 


±0.9 


14.7 


±0.9 


1.9 


0.24 


±0.06 


b 


980703 


BATSE 




0.966 


(8) 


14.15 


±0.1 


7.19 ± 


0.05 












0.11 


±0.01 


b 


990123 


BATSE/GRBM/Konus 


WFC 


1.600 


(9) 


19.81 i 


:0.03 


7.62 ± 


0.01 


35.6 


±0.2 


13.70 


±0.08 


1.8 


0.30 


±0.05 


b 


990506 


BATSE/GRBM/Konus 




1.3066 


(10) 


3.83 ± 


0.02 


1.66 ± 


0.01 












0.20 


±0.05 


n 


990510 


BATSE/GRBM/Konus 




1.619 


(11) 


2.54 ± 


0.03 


0.97 ± 


0.01 












0.20 


±0.01 


n 


990705 


GRBM 


WFC 


0.86 


(12) 


14.3 i 


:0.2 


7.7 ± 


0.1 


22.3 


±0.4 


12.0 


±0.2 


1.6 






b 


990712 


GRBM 


WFC 


0.433 


(13) 


4.1 ± 


0.2 


2.85 d 


:0.1 


4.8: 


±0.2 


3.35 


±0.1 


1.2 






n 


991208 


Konus 




0.7055 


(14) 


3.67 ± 


0.04 


2.15 ± 


0.02 
















n 


991216 


BATSE/GRBM/Konus 




1.02 


(15) 


3.80 + 


0.02 


1.88 ± 


0.01 












0.18 


±0.02 


n 


000131 


BATSE 




4.500 


(16) 


5.77 ± 


0.08 


1.05 ± 


0.01 












0.21 


±0.06 


n 


000210 


GRBM/Konus 


WFC 


0.846 


(17) 


2.4 ± 


0.2 


1.3 ± 


0.1 


5.35 


±0.3 


2.9: 


±0.2 


2.2 






n 


000214 


GRBM 


WFC 


0.47 


(18) 


2.5 ± 


0.4 


1.7 ± 


0.3 


6.8: 


±0.3 


4.65 


±0.2 


2.7 






n 


010222 


GRBM/Konus 


WFC 


1.477 


(19) 


3.68 + 


0.07 


1.48 ± 


0.03 


42.2 


±0.4 


17.0 


±0.2 


11.5 






n 


010921 


GRBM 




0.451 


(20) 


9.8 ± 


0.3 


6.75 h 


:0.2 
















b 


011121 


GRBM/Konus 


WFC 


0.362 


(21) 


10.0 d 


:0.3 


7.35 h 


:0.2 


18.9 


±0.2 


13.9: 


±0.15 


1.9 






b 


030329 


Konus 




0.1685 


(22) 


2.6 ± 


0.1 


2.19 ± 


0.08 
















n 



The energy range of BATSE data which was taken as reference instrument. For Konus and GRBM data the actual energy ranges are 50-200 

ke V and 40-700 keV re spect ively. 

(1) iBloom etail ll200ll); (2) iMeitzger et alj dl997h: (3) [Dioreovski etafl ll200ll): (4 ) [Kulkarni et al] ill 9981): (5 ) iBloom et alj d 19991): (6) 
iLamb, Castander & Reichartl dl999l); (7 ) iTinnev et al.l \l99& ~(8) Diorgovski et alj Jl998h: (9) iKulkarni et alj dl999h: (10 ) iBloom et alj 
d2003l): (11) iBeuermann et alj dl999h: (12) lAmati et all feoOOl); ( 13) IVreeswiik et alj feOOll) ; (14) iDodonov et alj Il999l) ; (15) 
Vreeswiik et alj dl999h: (16) lAndersen et all d200Ch; (17) iPiroet alj J2002h : (18) Untonerii et all d200Ch : (19) Ijhaetalj feOOlh : (20) 
Price et alj d200ai ; (21) lGarnavich et alJd2003l) : (22) lGreiner et al.ld2003h . 



Index <f has been estimated for an energy window width of two BATSE channels, 
(n) and (b) indicate narrow and broad width ACF class, respectively. 



gered LADs (usually the four closer to the line of sight). This 
combined data format was used when available, since the con- 
catenated pre-burst data allows a better estimation of the back- 
ground. In the case of GRB 970828 the DISCSC data are in- 
complete, and we used instead the 16-channel MER data type, 
binned up into 4 DISCSC-like energy channels. All BATSE 
bursts with known z were considered for study, excluding two 
cases were the data are incomplete or were no recorded (i.e., 
GRB 980326 and GRB 9806 1 3), resulting in a total of 1 1 cases. 



We also include the set of bursts selected in B04 that were 
observed by Konus, which is a G RB detector on board the 
Wind mission dAptekar et al.lll995l) . Light curves of its bursts 
are publicly available at 64 ms resolution in the 50-200 keV 
energy band. The collecting area of this experiment is about 
20 times smaller than the one of BATSE and consequently, in 
most cases, the signal is too weak for our temporal studies. 



However, in lB04l a comparative ACF analysis of the bursts ob- 
served by both instruments showed a good agreement for very 
bright bursts, and the selection criteria for Konus cases were set 
requiring peak count rates larger than 3000 counts s" 1 and the 
availability of post-burst data, resulting in the 5 Konus GRBs 
included in B04 burst sample. 

The BeppoSAX mission, that operated between the years 
1996-2002, had broad energy band capabilities thanks to the 
combined operation of seve ral instruments. The Gamma Ray 
Burst Monitor (GRBM; see iFrontera et al.lll997l) covered the 
40-700 keV energy range, roughly matching the range of 
Konus and BATSE (i.e., 55-320 keV using 2 + 3 channels). 
Note that since the ACF and the PDS are quadratic functions 
of the number of counts, and generally there are more counts 
at lower energies, the agreement for these temporal analysis 
functions will depend mainly on having a similar lower-end en- 



4 



L. Borgonovo et al.: Variability classes in GRBs 



[BAT5£>GFTB390123 




1/ 




Fig. 1. Examples of GRB light curves from the our sample of 
bursts with known redshifts, where the time t' is calculated at 
the source rest-frame. The light curves on the left (right) panels 
have broad (narrow) ACF widths. There are no obvious mor- 
phological differences between the two classes, both presenting 
cases of simple and complex structures. 



ergy limit. We made use of the standard high resolution 7.8125 
ms LCs, but we binned them up into a 62.5 ms time resolu- 
tion to improve the S/N ratio. This has a negligible effect on 
the measurement of the time scales that concern us here and 
it reduces some noise artifacts (i.e., the noise becomes more 
Poissonian). In addition, it approximately matches the standard 
64 ms BATSE temporal resolution for better comparison. The 
LCs were dead-t i me co rrected and background subtracted. In 
Borgonovo et alj ( 12005b a comparative ACF analysis was pre- 
sented, including 17 GRBS detected by the GRBM with known 
redshift from which 8 were also observed by BATSE. It was 
concluded there that although the measured dispersions of the 
local ACF widths were larger than in the BATSE case (~ 15% 
at half-maximum), the average sample values for each ACF 
width class were equal to within uncertainties. Therefore, at 
least in the context of the present temporal variability analysis, 
we can consider that these data have more intrinsic dispersion, 
but the introduced errors are mainly stochastic. 

The two Wide Field Cameras (WFCs) also on board 
Bepp oSAX covered the 2-26 keV energy range djager et al 



119971) . During their operation time they detected 53 GRBs in 
conjunction with the G RBM allowing the first broad band stud- 
ies of GRBs (see, e.g., Amati et al]|2002 ). Furthermore, a con- 
siderable fraction of the WFC bursts (~ 36%) were also de- 
tected by BATSE. The LCs were extracted with a time resolu- 
tion of 62.5 ms and discriminated in three energy channels (i.e., 
2-5, 5-10, and 10-26 keV). The energy intervals were chosen 
in order to have a similar amount of counts in each channel 
for a typical GRB. However, except for the brightest GRBs, 
the signal in each channel is too weak for the purposes of our 
temporal analysis and the LCs had to be integrated into a sin- 
gle energy channel. Here we will focus on the analysis of the 
13 GRBs for which we have redshift estimations. These bursts 
constitute a subset of our sample of GRBs with known redshift 
in the y energy band. Table [TJ lists all these bursts indicating 
in its first columns their name, the source instruments in y and 



X-ray bands when available, the estimated redshift z, and the 
corresponding reference. 

To further study the ACF width energy dependence we se- 
lected from the BATSE current catalog all long bursts (i.e., du- 
ration time T()o > 2 s) with a peak flux measured on the 1 .024 s 
timescale Fi s > 4 photons cirT 2 s~' in the 50-300 keV band 
that have available concatenated 64 ms data. This resulted in a 
sample of 188 bright bursts for which, in most cases, the red- 
shift is unknown. 

3. Methods 

For the autocorrelation function analysis we follow the same 
met hod presented in IB04L that was based on earlier works 



of iLink et alJ (119931) an d IF95L Here, we will summarise the 
method, and we refer to B04 for further details. Following the 
same notation, from a uniformly sampled count history with 
AT time resolution and N time bins, let m,- be the total observed 
counts at bin i. Also let bj be the corresponding background 
level and c, = m, - fr, the net counts. The discrete ACF as a 
function of the time lag t = kAT is 



N-l 



A(t = kAT) = J] 



, k = Q,...,N-l, (1) 



where 6 is the Kronecker function. Here the periodic boundary 
conditions (c, = c,- + #) are assumed. The normalisation constant 
Aq is defined as 



N-l 



(2) 



such that A(0) = 1 for k — 0. The term m, in Eq.|2]subtracts the 
contribution of the uncorrelated noise assuming that it follows 
the Poisson statistics. 

For practical reasons, the actual calculation of Eq. [TJ was 
done using a Fast Fourier Transform (FFT) routine. Denoting 
by Cf the Fourier transformed of the background subtracted 
light curve, then the definition of the power density spectrum 
(PDS) can be written as Pj = |C/| 2 . The noise contribution 
is subtracted from the PDS assuming Poisson statistics. Since 
the PDS and the ACF are Fourier pairs (Wiener-Khinchin the- 
orem), the latter is obtained by inverse transforming the first. 
Zero padding of the time series was used to avoid the artifacts 
produced by the periodic boundary condition. The normalisa- 
tion used for the ACF, that gives to its central maximum unity 
value, is equivalent (aside from a noise correction term) to the 
scaling of the LC by the square root of total power VPtot, where 
Aot = c r This normalisation is a natural choice for the ACF 
analysis and it makes the ACF of each burst independent (to 
first order) of its brightness. Note however that for our PDS 
analysis we found it more suitable to scale the LCs by their re- 
spective net count fluences T = c, (or equivalently dividing 
the PDS of the original LCs by T 2 ). Since the ao zero order 
Fourier coefficient is equal to the fluence T, our normalisation 
makes the PDS converge to unity towards the low frequencies, 
therefore in Fourier space it can be interpreted also as normalis- 
ing by the amount of power at the lowest frequency. In Sect. 14. 21 
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we will further discuss this choice over other possible normal- 
isations. 

The background estimations were done by fitting with up 
to a second order polynomial the pre- and post-burst data, that 
were judged by visual inspection to be inactive. This is critical 
for weak bursts and particularly aggravated in the Konus case 
since the publicly available LCs have almost no pre-burst data. 
This is the main reason why dim bursts were excluded in our 
sample selection in Sect. [2] 

In Ib04| it was found empirically, considering the analysis 
of the average ACFs of the narrow and broad width sets and 
their dispersion, that at half-maximum the separation between 
the two sets is most significant. Therefore, the ACF half-width 
at half-maximum w given by A(t = w) = 0.5 (here we deviate 
from the B04 notation for simplicity) was chosen as the mea- 
sure that best characterises each class. Since it has been shown 
that the average ACF decreases approximately following a 
stretch exponential dF95l) . the ACF width w was calculated fit- 
ting the logarithm of the ACF in the range 0.4 < A(r) < 0.6 
with a second degree polynomial. 

The uncertainties due to stochastic fluctuation were esti- 
mated using a Monte Carlo method. For any given LC a large 
number of realisations is generated assuming that the fluctu- 
ations on the gross number of counts in each bin m, follows 
a Poisson distribution. Then, the LCs are background sub- 
tracted and the ACF widths w are calculated following the same 
method as for the original data. However, when subtracting 
the Poisson noise contribution in the ACF (Eq. |2]l a factor of 
2 must be introduced in the corresponding noise term to take 
into account that the procedure to generate the synthetic LCs 
doubles the noise variance, since they are in a sense second or- 
der realisations. In general, if one considers a Poisson process 
with expected value p and iteratively assumes its counts as ex- 
pected counts of another Poisson process, after n iterations the 
expected value is still p, but the variance is np (with n — 2 in 
this case). Finally, the standard deviation of the obtained widths 
is used as estimator of the uncertainty. Nevertheless, the main 
source of uncertainty in dim bursts is the background estima- 
tion, that introduces a systematic error in the determination of 
the ACF width. Therefore, the reported uncertainties should be 
regarded as lower limits. Another likely source of systematic 
errors is the non uniform response of the detectors, since we do 
not deconvolve the count data. 

4. Analysis of the y-ray LCs 

Some LC examples from our sample of GRBs with known red- 
shift are shown in Fig.Q] No morphological differences are ev- 
ident between the narrow and broad classes by simple inspec- 
tion. However, standard linear analysis tools reveal clear vari- 
ability differences between the two, as we report in this section. 

4.1. ACFs 

With the inclusion of GRBM bursts we were able to expand 
the sample presented in IboH The added ACFs (shown in Fig. [2] 
in dashed lines) also follows a bimodal distribution when cor- 
rected for time dilation effects, reinforcing the previously found 




Fig. 2. a) Autocorrelation functions (ACFs) of 22 GRBs with 
known redshifts calculated in the observer's frame. Although 
easily identified in the rest-frame, the narrow (gray lines) and 
the broad [solid lines) classes overlap giving the impression of 
a unimodal distribution, b) Local ACFs where the cosmic time 
dilation effect has been corrected, being r' = r/(l + z). The 
newly added GRBM data (dashed lines) reinforce the bimodal 
pattern previously found by lB04l that used only the BATSE and 
Konus data (here both shown with solid lines). 



pattern. In Fig. [2^ the observed ACFs A(t) are shown first 
for comparison, while Fig. \2)p shows the rest-frame or local 
ACFs A(t'), where r' = t/(1 + z) is the time lag corrected for 
cosmic dilation. We obtained mean values for the rest-frame 



ACF half-width at half-maximum w 



(1.56 + 0.15) sand 



u" 



,(b) 

y 



= (7.42 ± 0.14) s, and also sample standard deviations 
o-y' = 0.6 s and cry = 0.5 s for the narrow and broad sub- 
se ts res pectively. The mean values are equal to those reported 
in IB04| within uncertainties. The gap between the two subsets 
(defined as the mean difference) represents a 4.9<x separation, 
smaller than the 7cr separation reported in lB04l mainly due to 
the increased dispersion in the broad class. However, thanks 
to a larger sample, we have slightly increased the significance 
of this separation. The probability p of a random occurrence of 
such gap was estimated using Monte Carlo methods. Assuming 
that there are no characteristic timescales and an underlying 
uniform probability distribution (the most favourable case) we 
estimated p < 4 x 10~ 7 . 
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Fig. 3. Power density spectra (PDS) of the subset of 15 bursts 
with narrow width ACFs. The frequency scale has been cor- 
rected for cosmic time dilation effects and the noise level has 
been subtracted assuming Poisson statistics. The upper panel 
shows all individual PDS where the frequency data have been 
equally binned in the logarithmic scale. The lower panel shows 
the estimated sample median {solid line) and the quartiles about 
the median {dashed lines) to indicate the dispersion. The decay 
phase is well modelled by a power-law {gray line) with index 
1.97 ± 0.04. Also shown, a Lorentzian function {dotted line) 
provides a fairly good fit over the whole frequency range. 

Hereafter our sample of bursts with known redshifts in the y 
band will be divided for analysis into two subsets of 15 narrow 
and 7 broad ACF width bursts respectively. In Table[T]we show 
the n ewly added ACFs widths together with those presented in 
B04l for completeness. 



4.2. PDS 

We calculate the PDS for our sample of GRBs with known red- 
shift, correcting the LCs for cosmic time dilation effects. We 
group bursts following the classes established in § !4.1l based on 
the ACF width. The upper panels of Figs. [3] and |4] show for the 
narrow and broad classes all individual PDS overlaid for com- 
parison. Given the chosen normalisation, all PDS must con- 
verge to unity at low frequencies (see Sect. [3j. The frequencies 
have been equally binned on the logarithmic scale. The proce- 
dure not only smooths out the PDS stochastic variations but en- 
able us to estimate an average PDS for each class. This is neces- 



sary because even if all the LCs had the same time duration, the 
redshift correction would make the Fourier frequencies differ 
for each burst (i.e., since the observed frequency / transforms 
to the rest-frame as /' = /(l + z))- The lower panels of Figs. [3] 
and|4]show the corresponding central values Pf, where we have 
chosen the use of the median over the mean as a more robust 
estimator of the expected or "underlying" spectrum. Although 
qualitatively similar results are obtained in both cases, the me- 
dian shows smaller fluctuations and it is less sensitive to the 
random exclusion of a few bursts from each class. The figures 
show also the quartile deviation {dashed lines) as a measure of 
the dispersion that was chosen to be consistent with the use of 
the median (i.e., the 25% and the 75% quartiles around the me- 
dian value). Comparing the two median PDS it is evident that 
they have remarkably different shapes. The narrow class PDS 
as expected shows more power at high frequencies and it is well 
described by a single power-law model Pj cc l/f a up to a low 
frequency cut-off at a; 0.1 Hz. In Fig. 03 we show a fit to the 
median PDS in the frequency range 0. 1 < f < 10 Hz for which 
we found a best-fit-parameter ay = 1.97 ± 0.04 (i.e., y band - 
narrow class). The obtained exponent is consistent with that of 
Brownian or red noise (i.e., with a = 2), which suggests the al- 
ternative use of a single Loren tzian model to describe the entire 
PDS as in shot noise models dBellilll992h . However, as shown 
also in Fig.[3]3, there seems to be a slight systematic deviation 
around the cut-off region. The broad class Pf on the other hand 
appears to have two distinct components, i.e., a broad low fre- 
quency component with a sharp break and a power-law compo- 
nent (Fig.|4j)). If we assume that the two components are inde- 
pendent then the single power-law should have an independent 
cut-off at low frequencies which cannot be determined but only 
constrained by our data. The low frequency component is well 
fitted by a stretched exponential P\{f) = exp[(-///i)'']. The 
best-fit-parameters for this model will depend to some extent 
on the cut-off frequency parameter of the power-law compo- 
nent and vice versa, however this will affect mainly the estima- 
tion of the exponent rj. Taking into account these uncertainties 
we found a characteristic frequency f\ = (0.025 ± 0.003) Hz 
and an index 77 = 1.3 ± 0.2. We considered other empirical 
models, concluding that functions with power-law asymptotic 
behaviour (e.g., a Lorentzian) do not decay fast enough to fit 
the data. For the high frequency component we found a best-fit 
power-law index or*' = 1.6 + 0.2 (i.e., y band - broad class), 
consistent with the Kolmogorov spectral value 5/3. However, 
note from Fig. [4^ that for different bursts the latter component 
appears to vary much more than the first. 

As we discussed in Sect. [T] the obtained mean PDS (as 
well as the median) may depend on the chosen normalisation 
since each burst will be weighted in a different way. If all LCs 
in the sample are realisations of the same stochastic process 
then the norm should simply improve the convergence to the 
mean value but should give identical results for a sufficiently 
large sample. Aside from the fluence norm, we tested several 
ways to scale the LCs, e.g., the peak flux used previo u sly by 
iBeloborodov et alj (1 19981 12OO0I) , the iMivamoto etail d 19921) 
normalisation (which expresses the PDS in fractional root- 
mean-square), and the root of the total power V^tot used for the 
ACF. We found in all cases qualitatively the same median PDS 
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Fig. 4. Power density spectra (PDS) of the subset of 7 bursts 
showing broad width autocorrelation functions (ACFs). As in 
Fig.[3]the individual PDS and the sample median are shown for 
comparison. The median PDS is significantly different to that 
of the narrow class and two components appear to be present. 
The low frequency component has a sharp cut-off and is well 
modelled by a stretched exponential {gray line). The high fre- 
quency component is well described by a power-law with index 
1.6 + 0.2 consistent with a 5/3 Kolmogorov spectral index. 



for each class but significantly larger spreads. Obviously, this 
does not prove that each class sample is "uniform" since it is a 
necessary but not sufficient condition. In this respect, it is worth 
mentioning that when combined into a single sample and with- 
out correcting for redshifts effects we obtained a median PDS 
(shown in Fig- [5j that is well fitted in the 0.01-1 Hz frequency 
range by a single power-law with index 1.73 ± 0.07, consistent 
within uncertai nties with the Kolmogorov spectrum and in full 
agreement with lBeloborodov et al.l d 1998b . Therefore, our sam- 
ple of GRBs with known z is not different, at least in this regard, 
to the general sample of bright BATSE GRBs used in previous 
works. Note also that the overall dispersion about the median 
is significantly larger than in Figs. [3] and |4] where, despite the 
smaller sample sizes, redshift corrections and ACF classes have 
been taken into account. 
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Fig. 5. Median power density spectra of our whole sample 
of 22 GRBs with known redshift calculated in the observer's 
frame. The spectrum shows a power-law behaviour approxi- 
mately within the 0.01-1 Hz frequency range and a best fit (dot- 
ted line) gives an index 1 .7 3±0.07, consistent with the results of 
Beloborodov et alj d 1 998b derived without redshift corrections. 
Note that the quartile dispersion (dashed lines) is considerably 
larger than those found in Figs. [3] and [4] 



into the X-ray energy range. These bursts represent a subset of 
the previously used sample in the y-ray band (Sect.[4]i. We find 
also at these energies significant differences between the two 
ACF classes although the results inevitable have larger associ- 
ated uncertainties. 



5.1. ACFs 

Although the ACF widths are substantially broader at these 
lower energies (2-26 keV), the two ACF classes are still easily 
distinguishable as shown in Fig. [6] However, the GRB 010222 
classified as having a narrow ACF in the y band (gray lines) has 
broaded to such extent that it falls into the broad width range, 
although its ACF decays more slowly than any of the broad 
width cases. As shown in TableQ]where we list all the observed 
and the local ACF widths, this is the only outlier in the sam- 
ple. In column 10 we list as a measure of the broadening the 
ACF width ratio for the y and X-ray bands. While the average 
width ratio is (wx/w r ) 2, the broadening in the GRB 010222 
case (wx/wy 1 1.5) is much larger than any other burst in the 
sample and it is most likely caused by a systematic error, as 
we will discuss later in Sect. [7] For this reason this case will be 
kept separate from the sample for the rest of our WFC data tem- 
poral analysis. Consequently, the obtained average local width 



ACF are w 



,(n) 



(3.1 ± 0.5) s and w 



,(b) 



(13.6 ± 0.6) s, with 



sample standard deviations (relative dispersions) cr 



.(ID 



1.9 s 



.(b) 



(60%) and a x 
respectively. 



1.6 s (12%) for the narrow and broad subsets 



5. Analysis of the X-ray LCs 

In this section we use the WFC sample of 13 GRBs with known 
redshifts (see TableHJ to extend the previous temporal analysis 



5.2. PDS 

Once again using our sample of WFC bursts with known z, we 
estimate the local median PDS for each ACF width class fol- 
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lowing the same methods used in Sect. 14.21 We consider sepa- 
rately the special case GRB 010222 due to the unique broaden- 
ing of its ACF. Figure [7] shows our results for the narrow class. 
As expected from the ACF analysis the break appears now at 
approximately half the corresponding frequency at y energies 
(Fig. [3j, but most noticeably the behaviour towards the high 
frequencies now shows a much faster decay. A best-fit for fre- 
quencies /' > 0.08 Hz using a Pf oc 1 jf a model gives an expo- 
nent af = 3.0 ± 0.2, a significantly steeper power-law than in 
the y band. The PDS of GRB 010222 (gray line) is also shown 
in Fig. [7] for comparison. Although with much less power at 
low frequencies than the other bursts in this class, at high fre- 
quencies (/' > 0.2 Hz) it follows a very similar asymptotic 
behaviour. 

Our estimation of the median PDS P / for the broad class 
is shown in Fig. [8] In spite of the large uncertainties associ- 
ated with such a small sample, the two components found in 
Fig. |4] are recognisable. The low frequency component is not 
as prominent and broad as before, and consequently the data 
are not enough to properly constrain the fit parameters of a 
stretched exponential as it was previously done. However, we 
estimated that the high frequency component follows a power- 
law with index — 1.7 ± 0.2 consistent within uncertainties 



with the Kolmogorov index found in Sect. 14.21 for the y-ray 
LCs. When compared with the broad class Pf in Fig. [8] the 
GRB 010222 spectrum matches neither the Pf general shape 
nor any of its components. Power-law fits using the indices and 
the frequency ranges shown in Figs. [7] and [8] for the respective 
high-frequency components resulted in squared residual totals 
per degrees of freedom x 2 l v = 16.2/11 andx 2 /v = 179. /17 
for the narrow and broad PDS respectively. Based on these fits 
we conclude that the first model is acceptable while the sec- 
ond model can be rejected at a very high confidence level. 
Uncertainties were estimated using synthetic shot noise and 
calculating standard deviations on a sample of PDS follow- 
ing the same procedures and frequency binnings applied to the 
burst data. 

6. Comparing ACF classes 

In order to understand the origin of the two classes found 
based on our temporal analysis we looked for any additional 
GRB characteristic or physical parameter that might differ be- 
tween them. As point out in Sect. |4j the visual inspection 
of the LCs reveals no trivial morphological differences be- 
tween the classes (see Fig. Q}. Both classes present cases with 
simple and complex structure, showing from a few smooth 
pulses to many heavily overlapped sharp pulses. To cover 
a broad range of physical parameters, we made use of the 
databas e of GRB redshifts and o ther burst parameters com- 
piled in Fried man & Bloom ( 20051) and references therein, that 
contains much additional information about the 22 bursts in 
Table[U although not all parameter estimations are available for 
every burst. In particular, we considered the following parame- 
ters (primed quantities are calculated at the source rest-frame): 
the peak energy E' pk , the isotropic equivalent gamma-ray en- 
ergy Ei so , the duration time TL, the time of the observed break 
in the afterglow light curves 77, and the redshift z. In all cases 




Fig. 6. Autocorrelation functions (ACFs) calculated at the 
source rest-frame for a sample of 13 GRBs detected by WFC 
in the 2-28 keV energy band. In spite of the broadening of the 
ACF at lower energies, the bimodal distribution found in the y- 
ray band is clearly seen, although with larger relative dispersion 
for each class. However, one burst (GRB 010222) previously 
classified as belonging to the narrow ACF width class (gray 
lines) appears wider than the broad ACF width class (solid 
lines). This behaviour is most likely not intrinsic but due to 
a systematic error (see discussion in Sect. 13. 
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Fig. 7. Median power density spectra (PDS) from the subset 
of WFC bursts with known redshifts belonging to the narrow 
width ACFs class, excluding the outlier case GRB 010222 that 
is shown along side (gray line). The decay at high frequencies 
(/' > 0.1 Hz) approximately follows a power-law, but with a 
significantly steeper slope than the previously found in the y- 
ray energy band (3.0+0.2). GRB 010222 shows similar asymp- 
totic behaviour at high frequencies. 



we compared the distributions within the narrow and broad 
ACF classes using the standar d Kolmogorov-Smirnov (K-S) 
test (see, e.g. JPress et al.l[l992l) . which is particularly sensitive 
to median deviations, but we found no significant difference 
for the cumulative distributions of any of those parameters. No 
significant R correlation coefficients with w' were found for £7 k 
and E[ so (logarithms taken in all cases), although marginal re- 
sults were obtained for the temporal parameters TL (R = 0.44 
with significance p < 0.04) and (R - 0.46 with significance 
p < 0.05 for only 18 bursts with determined values). Note 
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ditional redshift corrections to the accounted cosmic dilation 
were studied using the transformation 



0.05 0.1 
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Fig. 8. Median power density spectra (PDS) from the subset 
of 4 WFC bursts with known redshifts belonging to the broad 
width ACFs class. Despite the small sample size, the two broad 
PDS components identified at y-ray energies are recognisable. 
The index of the power-law component (dotted line) is consis- 
tent with the Kolmogorov value as previously found in Fig. [4] 
given a best-fit-parameter 1.7 + 0.2. 



that a correlation between w and Tqq has been well established 
(R = 0.58), and although in the observer-frame is partially due 
to the cosmic temporal dilations, it should b e found also in 
the rest-frame (IBorgonovo & Biornssonl 120061) . However, our 
present sample is too small to determine it at a high confidence 
level. More surprising is the finding of an indication of some 
correlation between a prompt (w') and an afterglow (7V) pa- 
rameter, but this would need better statistics to be confirmed. 

Other subdivisions of the long/soft GRB class have been 
suggested in the literature but we found them to be unrelated 
to with the one under consideration. Evaluation of the distri- 
butions of the observed and E v \ for the broad ACF class 
excludes an association between these bursts and the "inter- 
medi ate" duration GRB gro up reported by iHorvafhl (1998) 



and iMukherjee et al 



(1998). Likewise, we found no corre- 
spondence between the ACF bimodality and the bimodality 
of the afterg l ow op tical luminosity dis t ributio n reported by 



oi tne aiterg l ow op tical luminosity a istnbutioi 
Nardini et all ( 120061) and lLiang & Zhangl d2006l) . 



7. ACF energy dependence 

In our analysis we have only corrected the LCs for cosmic time 
dilation. However, the data are collected over a finite energy 
band and therefore additional corrections might be needed to 
account for the shift in energy, since in principle we should 
compare LCs emitted in the s ame energy band. As LC pulses 
are sharper at higher energies (INorris et al]l 19961) . correspond- 
ingly the ACF is narrower. For a sample of 45 bright bursts, 
lF95l found that the width of the mean ACF depends on the 
energy E as w(E) oc E~ 0A3 . Since for large redshifts the in- 
strument will see photons emitted at higher mean energies, the 
anti-correlation of the ACF width with the energy should par- 
tially counteract the time dilation effect. In lB04l possible ad- 



w/(l +zY 



(3) 



for the observed widths, where a represents a smaller than unity 
correction. The relative dispersion of each width class was cal- 
culated as a function of the parameter a looking for minimum 
values (see Fig. 6 in |B04|) since, under the assumption that the 
W intrinsic distribution does not depend on z (i.e., evolution 
effects are neglected), redshift dependencies should increase 
the observed dispersion. We obtained minimal dispersions at 
a (n) = -0.35 ± 0.2 and a (b> = 0.04 + 0.04 for the narrow and 

min min 

broad classes respectively, equal within uncertainties to previ- 
ous results. The uncertainties w ere estimated num erically using 
the bootstrap method (see, e.g., Press et aTT l992). If the disper- 
sion is calculated over the whole sample (i.e., without separat- 
ing the two width classes) no well-defined minimum is found, 
with an approximately constant value within the free param- 
eter range |a| < 1, which further supports the separation into 
two classes. However, while the narrow class minimum agrees 
within errors with the expected effect from the energy shift, the 
broad class dispersion does not improve much with any addi- 
tional redshift correction. Actually, in both cases a perfect can- 
cellation (i.e., amin = 0) cannot be ruled out, although it would 
appear coincidental that all other possible redshift dependen- 
cies would exactly compensate the energy shift. Another pos- 
sibility is that for the broad class the ACF dependence with 
energy is weaker. At face value the low dispersion already sug- 
gests this and in this case a nearly null a correction would be 
much more likely. In what follows we will show evidence that 
supports the latter alternative. 

Using 4 energy channel data from a sample of 188 bright 
BATSE bursts we analyse the ACF width dependence with en- 
ergy for each individual case using a single power-law model 
w(E) oc E~^ (the channels are approximately equally spaced 
on a logarithmic scale). We found that in most cases the rela- 
tion w(E) in individual bursts is not well fitted by a power-law 
(less than 33% of the bursts satisfied for the sum of the reduced 
residuals x% < 2), but for most cases it can be considered a 
good approximation within the BATSE energy range with typ- 
ical total relative errors < 10%. The obtained distribution of £ 
indices is shown in Fig. [9] Since we found for our fits poorer 
figures-of-merit (i.e., xl) f° r larger £ values, the largest values 
are related to the cases that strongly deviate from the power-law 
model for w(E) and are therefore ill determined. For this same 
sample, a power-law fit of the energy dependence of the width 
of th e mean ACF gives an index -0.40 ±0.03 in agreement with 



F95l However, note in Fig.[9]fhat for such a long-tailed distri- 



bution most individual bursts typically have a weaker energy 
dependence and therefore the median fi c h = 0.29 is a more 
representative value of individual indices. 

In addition, since we are particularly interested in quan- 
tifying the effect of the energy dependence in our ACF anal- 
ysis where the observation energy window has the width of 
two BATSE channels, we studied the ACF width energy de- 
pendence w(E) combining the 4-channel data into three over- 
lapped 2-channel width energy bands. The distribution of the 
thus obtained indices ^2ch is shown in Fig. [9] (dashed lines). 
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Fig. 9. Distribution of the index £ in the power-law relation 
w(E) oc used to model the energy dependence of the ACF 
width, obtained using 4-channel concatenated data (solid line) 
for a sample of 188 bright BATSE bursts. The distribution has 
an estimated median £i c h = 0.29 and it shows an approximately 
exponential decay towards larger £ values. For the purposes of 
our ACF analysis that uses the two central BATSE channels, we 
modelled also the w(E) relation using a 2-channel energy win- 
dow for each data point. As expected, the ACF width shows a 
weaker energy dependence and consequently lower typical val- 
ues for the £ index. The thus obtained distribution [dashed line) 
is narrower and with a significantly lower median %2ch = 0.21. 



The median §2ch = 0.21 is significantly smaller, as expected 
since less correction will be needed the wider the energy win- 
dow. If we consider the median value of the ACF broadening 
wx/wy for the sample in Table [TJ (excluding GRB 010222) we 
derive a median index £ « 0.20, in good agreement with the 2- 
channel energy window distribution in Fig. [9] either taking the 
lower edge of each energy window (following |F95l) or their ge- 
ometric means. In Table Q] we list our estimations of the index 
£ 2c h for the 1 1 BATSE GRBs with known z. Both ACF width 
classes show typical values that seem consistent with the gen- 
eral distribution in Fig. [9] but unfortunately the small sample 
size only allowed us to distinguish very substantial differences. 

Since we do not have redshifts for most of the BATSE 
GRBs, the estimated distribution of index £ was derived mix- 
ing both ACF width classes; nevertheless it is probably fairly 
representative of the narrow class that dominates the sample. 
In principle, we do not know how representative that general 
distribution is of the broad class. For our sample of 1 1 BATSE 
GRBs with known z the fraction of broad cases is ft, = 4/11 
(/b = 7/22 for our whole sample). However, since the after- 
glow localisations needed for the redshift determinations al- 
ways involved another instruments with different trigger re- 
sponses, that fraction is most likely not representative of the 
general BATSE catalog. An analysis of the ACF widths of the 
188 brightest long BATSE bursts sample shows that the dis- 
tribution is well described by a log-normal distribution with 
a mean 10* log1 ^ = 3.0 s and a standard deviation represent- 
ing a factor of 3.1, with no significant hint of an underlying 
bimodality (see also iBorgonovo & Bjornsson 2006). Since all 
1 1 BATSE GRB with known redshift have w > 2.5 s almost 



half of the BATSE w distribution without redshifts seems trun- 
cated. Using the K-S test that compares cumulative distribu- 
tions we calculate a chance probability p < 0.02 that they are 
drawn from the same distribution. We conclude that the sample 
with known z is significantly biased towards large ACF widths, 
mainly affecting the representation of the narrow class. Taking 
into account this truncation and the fact that for a burst to be- 
long to the broad class its observed ACF width must be w > 7 s 
independently of the redshift, we estimate that the unbiased 
BATSE fraction of broad cases should lie approximately be- 
tween 0.08 < fb ^ 0.22 to be consistent with the observations. 

Having characterised the energy dependence of the 
timescale w with the index we investigated possible corre- 
lations with other GRB parameters, but bearing in mind that 
the largest £ values are most likely meaningless. We tested a 
broad range of temporal and spectral parameters (e.g., duration 
time 7*90, emission time, ACF width w, time lag between en- 
ergy channels, vF v peak energy E p ^, low and high spectral in- 
dices a and &) following definitions and methods described in 
B org onovo & Bjornsson (2006) for a multivariate corre ation 



analysis. In particular, we used a measure of the overall spectral 
evolution of a burst introduced there, where the LC is divided in 
two fluence-halves and the peak energies of the corresponding 
integrated spectra are estimated and compared. In that way a ra- 
tio of peak energies is defined as fiE^ = £pk (1) / £pk <2 \ so that 
■REpk > 1 implies an overall hard to soft evolution. Although 
we found no significant global correlations with the index 
there are indications of a different behaviour at low £ values in 
two of the studied parameters, i.e., w and HE^. 

In Fig. [TO] we show scatter plots of the ACF width w and 
the ratio of peak energies 7?£pk versus the index The £ ob- 
served range was divided into four bins with equal number 
of data points, and the corresponding geometric means (more 
appropriate for log-normal distributions) and uncertainties are 
shown. Due to the large spread the variables seem uncorrected 
based on their correlation coefficients R, nevertheless the bin- 
ning reveals a trend where lower £ values correspond on av- 
erage to broader ACF widths w and lower values of the < R£ p k 
ratio. In addition, a weak but significant anti-correlation exists 
between w and flE p k with coefficient R = -0.25 taken their 
logarithms. For £ < 0.3 there is noticeably less dispersion in 
log 'REpk (i.e., less relative dispersion in < RE p k) and a clear pos- 
itive correlation that comprises half of the sample (with coef- 
ficient R = 0.45). Comparing a linear fit to those data to the 
no correlation null hypothesis gives for the F-test (the standard 
test for a ratio of x 2 values) a level of rejection p < 0.001. 
However, since we do not have an objective criterion to trun- 
cate the index f distribution, the real significance of this found 
correlation is hard to evaluate; therefore we approach the prob- 
lem using Monte Carlo methods. We draw random samples for 
both parameters using the bootstrap method so that they are 
uncorrelated. Then sorting the data based on their £ values, we 
look for the ^-truncate that would give the maximum chance 
correlation with HE^. Assuming conservatively that the corre- 
lation has to comprise at least 1/3 of the sample and \R\ > 0.3 
then the trial is judged as successful. In this way we found a 
probability of p < 0.008 of finding a similar or stronger corre- 
lation by chance. Although the energy dependence of the ACF 
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is clearly a consequence of the spectral variability, the spectral 
evolution of a burst is as complex a s its LC, shifting from har d 
to soft over each individual pulse (IBorgonovo & Rvdell200ll) . 
For this reason it is only under certain limited conditions that 
we are able to actually observe a correlation between £ and 
fiEpk. In summary, from the analysis of Fig. \W\ there are clear 
indications that bursts with broad ACF on average show weaker 
ACF energy dependence and less spectral evolution. 

Finally, we attempted to extend the analysis of the ACF 
energy dependence to X-ray energies with the use of 3-channel 
data from the WFC bursts which were simultaneously observed 
by BATSE. However, only the signal of the very brightest of the 
WFC bursts were found suitable for this purpose and the resul- 
tant sample was considered too small for a systematic broad- 
band analysis. Furthermore, although for most of these bursts 
w(E) shows a smooth, approximately power-law behaviour, we 
observed in a significant number of bursts a function disconti- 
nuity or "shift" between the WFC and BATSE data, a system- 
atic error most likely due to the background estimation. The 
observed fraction of cases that present this artifact (< 20%) 
is consistent with the occurrence of one such case within the 
1 3 WFC bursts with known z in Table Q] Therefore, we con- 
clude that the unusual ACF broadening of the outlier case 
GRB 010222 discussed in Sect. 15. II is most likely an artifact. 
The PDS analysis also supports this since an incorrect back- 
ground subtraction would mainly affect the power at low fre- 
quencies, and we have shown in Fig. [7] that GRB 010222 de- 
viates from the narrow class typical behaviour only at those 
frequencies. 

8. Discussion 



0.2 Q.i . ':' :■. £ 



1.2 



The reported ACF bimodality in B04 was based on the analy- 
sis of a small sample of GRBs with available z. Nevertheless, 
the gap between the two ACF classes was found to be highly 
significant (with p < 6 x 10~ 7 ). Low-number statistics can of- 
ten suggest nonexistent features and should be regarded with 
caution. As always, one should consider the possibility of a se- 
lection effect not accounted for in the probability estimation. 
However, while it is easy to conceive of upper and lower biases 
in the observed ACF width distribution (e.g., due to trigger cri- 
teria), it is harder to think of a reason why this would happen 
for middle range w' values. We have shown in Sect|7]that, when 
compared to a complete, flux limited large sample of BATSE 
bursts, our sample with known z is likely biased against very 
narrow ACFs. The addition of GRBM data with their larger in- 
trinsic dispersion only slightly increased the statistical signif- 
icance of the found ACF bimodality. Nevertheless, with three 
sets of y-ray data showing the same bimodal pattern (Fig. |2]i 
we can safely rule out an instrumental effect as the cause of it. 

It is clear from the PDS analysis presented in Sects. H] and 
[5]that narrow and broad ACF bursts have very different overall 
variability, with the latter showing a prominent low-frequency 
PDS component with a very sharp cut-off, which explains in 
part the associated low dispersion observed in this ACF width 
class, despite the typical fluctuations in individual burst PDS 
and the observed dispersion in the high-frequency power-law 
component. Another possible difference between the two ACF 




Fig. 10. Correlation analysis for the index £ from the relation 
w(E) cc E~s describing the energy dependence of the ACF 
width based on a sample of 188 BATSE bursts. Although we 
found no overall correlation between £ and other commonly 
used burst parameters (see Sect. [7]), at small £ < 0.3 values 
two parameters show a significantly different behaviour com- 
prising ~ 50% of the sample. At the lower corners we show 
average values of the estimated uncertainties representative of 
each data halves. Large £ values have much larger uncertainties 
and they are usually associated with poor fits. The £ axis has 
been divided into four bins with equal number of data points 
for averaging (solid squares). Upper panel: Typical values of 
the ACF width w are higher for lower £ indices. Lower panel: 
The ratio of peak energies 7?£pk, that is a measure of the burst 
energy spectral evolution, shows less dispersion and a positive 
correlation at low ^ indices. 



classes is hinted at by the trends shown in Fig.QI)] which could 
be interpreted as an indication that broad ACF bursts show less 
spectral evolution than the narrow ones. This view is indirectly 
supported by the comparison of the respective median PDS Pf 
at different energies estimated in Sects. 14.21 and [5721 The nar- 

(ft) 

row class exhibits a single power-law PDS with index at w 2 
at high energies and a steeper a£ ~ 3 low energies. In shot 
noise models where LCs are generated by random pulses, an 
index 2 commonly appears when the pulses have sharp rising 
phase, like the typical FRED-like pulses in many GRBs. Bu t 
due to the hard-to-soft evolution dBorgonovo & Rvdel l200lh 
the same pulses at X-ray energies look smoother with slower 
rise phases giving steeper indices. On the other hand, the high- 
frequency power-law component in the broad class PDS shows 
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approximately the same index in both energy bands (close to 
the 5/3 Kolmogorov index). Visual inspection in the y and X- 
ray bands of the brig ht broad ACF case GRB 990123 (see, 
e.g., Corsi et al . 2005), which shows well-separated pulses, in- 
dicates that the broadening of the ACF at low energies is not 
as much due to a corresponding smoothing of individual short 
pulses as to a decrease in their relative importance with respect 
to a long smooth LC component that seems to lie underneath 
them. If this case were representative of the whole class then 
that would explain in terms of LC components the observed Pf 
behaviour at different energies. 

One remarkable feature of the found ACF bimodality is the 
low relative dispersion of the broad class (~ 6% in the y band) 
around the w'*' = (7.42 ± 0.14) s timescale. If confirmed such 
a temporal feature would have very important implications, 
both for our physical understanding of GRBs and because of 
its potential use in cosmological studies. A crucial issue then 
is to properly account for redshift dependencies. The reported 
ACF widths were derived taking into account only the cosmic 
time dilation, since other effects are of a lesser order. The latter 
were quantified by means of the free parameter a included in 
Eq. [3] width transformation. We found well-defined minima at 
a (n) = -0.35 ± 0.2 and = 0.04 + 0.04 for the narrow and 

mm mm ■ ■ 

broad ACF class respectively (in agreement with lB04l) . Since 
|F95| reported a considerable energy dependence for the aver- 
age ACF (w(E) cc E~ 0A3 ) then the cosmic energy shift would 
be able to account within errors for the found a m i n in the nar- 
row class. However, for the broad class either we would need 
to assume that there are other dependencies at work that coin- 
cidentally cancel out almost exactly the energy shift effect or 
that simply the ACF energy dependence is weaker in this case. 
From the analysis of this alternative presented in Sect. [7] we 
conclude that in the context of our ACF study the power-law 
energy dependence is better characterised by a median index 
ftch = 0.21, still consistent with the narrow class within errors 
and a considerably smaller effect to account for. We evaluated 
by Monte Carlo simulations that the estimated w' dispersion 
(for a broad class sample of only 7 bursts) is smaller but not 
significantly different to the one expected assuming the %2ch 
distribution found in Fig. [9] Moreover, we found indications 
that the broad ACF class may have lower typical £ indices than 
the general sample. This would be an expected consequence of 
the weaker spectral evolution already mentioned above when 
comparing the two classes. 



Based on the analysis of WFC X-ray data, IVetere et al 



(2006) reported that in a fraction of bursts the existence of 
two LC components is apparent. They selected 10 bursts by 
visual inspection and found for their slow/fast (low/high fre- 
quency variability) LC decomposition that the slow compo- 
nents were in general softer than the corresponding fast ones, 
suggesting different underlying emission processes. There are 
only 3 GRBs that they claim to have slow components among 
our sample with known z, but they all show broad ACFs in the 
source rest-frame (including the outlier GRB 010222), which 
suggests a connection to our own results that will need further 
analysis to be established. T he discussion abou t the nature of 
the slow/fast components in IVetere et al.1 (§006) would apply 
also to the low/high frequency components found in broad ACF 



bursts. Wit hin the framework of the internal shock model de- 
veloped by M eszaros & Rees ( 2000l) . there would be two major 
radiating regions in the relativistic outflow and the PDS com- 
ponents would arise from the different expected variability of 
the emission, i.e., from a thermal photospheric component and 
from the non-thermal flux of optically thin dissipative regions 
above the photospheric radius. In this context, ACF class dif- 
ferences could be due to the conditions cont rolling the relativ e 
strength of those two emission components ( Ryd e et al . 2006). 



Acknowledgements. We thank S. Larsson and C.-I. Bjornsson for use- 
ful comments and careful reading of the manuscript. We are grate- 
ful to the anonymous referee for constructive criticism. This research 
has made use of BATSE and Konus data obtained from the High 
Energy Astrophysics Science Archive Research Center (HEASARC), 
provided by NASA's Goddard Space Flight Center. 

References 

Amati, L., et al. 2000, Science, 290, 953 

Amati, L., et al. 2002, A&A, 390, 81 

Andersen, M. I. et al. 2000, A&A, 364, L54 

Antonelli, L. A., et al. 2000, ApJ, 545, L39 

Aptekar, R. L. et al. 1995, Space Science Reviews, 71, 265 

Belli, B. M. 1992, ApJ, 393, 266 

Beloborodov, A. M., Stern, B. E., & Svensson, R. 1998, ApJ, 
508, L25 

Beloborodov, A. M., Stern, B. E., & Svensson, R. 2000, ApJ, 
535, 158 

Beuermann, K., et al. 1999, A&A, 352, L26 
Bloom, J. S., et al. 1999, Nature, 401, 453 
Bloom, J. S., Djorgovski, S. G, & Kulkarni, S. R. 2001, ApJ, 
554, 678 

Bloom, J. S., Berger, E., Kulkarni, S. R., Djorgovski, S. G, & 

Frail, D. A. 2003, AJ, 125, 999 
Borgonovo, L. 2004, A&A, 418, 487 (B04) 
Borgonovo, L., & Ryde, F. 2001, ApJ, 548, 770 
Borgonovo, L., & Bjornsson, C.-I. 2006, ApJ, 652, 1423 
Borgonovo, L., Frontera, F, Guidorzi, C, Montanari, E., & 

Soffitta, P. 2005, II Nuovo Cimento, 28 C, 275 
Corsi, A., et al. 2005, A&A, 438, 829 
Costa, E., et al. 1997, Nature, 387, 783 
Djorgovski, S. G, et al. 1998, ApJ, 508, L17 
Djorgovski, S. G, Frail, D. A., Kulkarni, S. R., Bloom, J. S., 

Odewahn, S. C, & Diercks, A. 2001, ApJ, 562, 654 
Dodonov, S. N., Afanasiev, V. L., Sokolov, V. V, Moiseev, 

A. V, & Castro-Tirado, A. J. 1999, GCN Report 475 
Fenimore, E. E., int't Zand, J. J. M., Norris, J. R, Bonnell, J. T, 

& Nemiroff, R. J. 1995, ApJ, 448, L101 (F95) 
Fishman, G. J., et al. 1989, in Proc. of the GRO Science 

Workshop, ed. W. N. Johnson, 2 
Friedman, A. S., & Bloom, J. S. 2005, ApJ, 627, 1 
Frontera, F, Costa, E., dal Fiume, D., Feroci, M., Nicastro, L., 

Orlandini, M., Palazzi, E., & Zavattini, G. 1997, A&AS, 122, 

357 

Garnavich, P. M. et al. 2003, ApJ, 582, 924 
Greiner, J., et al. 2003, GCN Report 2020 
Horvath, I. 1998, ApJ, 508, 757 
Jager, R., et al. 1997, A&AS, 125, 557 



L. Borgonovo et al.: Variability classes in GRBs 



13 



Jha, S. et al. 2001, ApJ, 554, L155 

Kolmogorov, A. 1941, Dokl. Akad. Nauk SSSR, 31, 538 

Koshut, T. M., Kouveliotou, C, Paciesas, W. S., van Paradijs, 

J., Pendleton, G. N., Briggs, M. S., Fishman, G. J., & 

Meegan, C. A. 1995, ApJ, 452, 145 
Kulkarni, S. R., et al. 1998, Nature, 393, 35 
Kulkarni, S. R., et al. 1999, Nature, 398, 389 
Lamb, D. Q., Castander, F. J., & Reichart, D. E. 1999, A&AS, 

138,479 

Lazzati, D. 2002, MNRAS, 337, 1426 
Liang, E., & Zhang, B. 2006, ApJ, 638, L67 
Link, B., Epstein, R. I., & Priedhorsky, W. C. 1993, ApJ, 408, 
L81 

Meszaros, P., & Rees, M. J. 2000, ApJ, 530, 292 

Metzger, M. R., et al. 1997, Nature, 387, 878 

Miyamoto, S., Kitamoto, S., Iga, S., Negoro, H., & Terada, K. 

1992, ApJ, 391, L21 
Mukherjee, S., Feigelson, E. D., Jogesh Babu, G, Murtagh, F, 

Fraley, C, & Raftery, A. 1998, ApJ, 508, 314 
Nardini, M., Ghisellini, G., Ghirlanda, G., Tavecchio, F, 

Firmani, C, & Lazzati, D. 2006, A&A, 451, 821 
Nakar, E., & Piran, T. 2002, MNRAS, 331, 40 
Norris, J. P., Nemiroff, R. J., Bonnell, J. T., Scargle, J. D., 

Kouveliotou, C, Paciesas, W. S., Meegan, C. A., & Fishman, 

G. J. 1996, ApJ, 459, 393 
Piro, L. et al. 2002, ApJ, 577, 680 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. 

P. 1992, Numerical Recipes in Fortran (2d ed.; Cambridge: 

Cambridge Univ. Press ) 
Price, P. A., et al. 2002, ApJ, 571, L121 

Reichart, D. E., Lamb, D. Q., Fenimore, E. E., Ramirez-Ruiz, 
E., Cline, T. L, & Hurley, K. 2001, ApJ, 552, 57 

Ryde, F, Bjornsson, C.-L, Kaneko, Y., Meszaros, P., Preece, 
R., & Battelino, M. 2006, ApJ, 652, 1400 

Tinney, C, et al. 1998, IAU Circular 6896 

Vetere, L., Massaro, E., Costa, E., Soffitta, P., & Ventura, G. 
2006, A&A, 447, 499 

Vreeswijk, P. M., et al. 1999, GCN Report 496 

Vreeswijk, P. M., et al. 2001, ApJ, 546, 672 



